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The motivation and the challenge in applying the renormalization group for 
Q-i! systems with several scaling regimes is briefly outlined. The four dimensional 

4 model serves as an example where a nontrivial low energy scaling regime 
is identified in the vicinity of the spinodal instability region. It is pointed 
out that the effective theory defined in the vicinity of the spinodal instability 
offers an amplification mechanism, a precursor of the condensation that can 



X 



be used to explore nonuniversal forces at high energies. 



I. INTRODUCTION 



The idea of the renormalization group |I[, is to view the different interactions in 
a hierarchical manner, by building up the more complex systems from their "elementary" 
constituents. The renormalization group should describe the manner the elementary parti- 
cles rearrange themselves in forming the composite particles captured by the detectors with 
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given resolution, the complex structure of ordinary matter and finally the transition from 
the micro to the macroscopic physics. Due to the obvious technical difficulties of such an 
ambitious project the realization of this idea is severely restricted. What is usually achieved 
by the help of different analytical approximation methods is to analyze the dependence on 
the observational scale in a scaling regime where the evolution equations are linearizable or 
at least perturbative. The result is the renormalized trajectory, the scale dependence of the 
coupling constants in a regime that is dominated by a given interaction. One may call this a 
local analysis of the renormalization group flow, performed in the individual scaling regimes 
that usually but not necessarily agree with the vicinities of the fixed points. A number of 
important results have been derived in this manner. 

The real challenge, we believe, is to describe the transmutation of one set of scaling laws 
into another one as the scale of the observation is changed. This requires the construc- 
tion of the renormalized trajectory connecting different scaling regimes. Such a manifestly 
nonperturbative phenomenon will be studied in this paper in the case of a simple model 
with two scaling regimes, the single component 4 scalar field theory in the spontaneously 
broken phase. This model supports an asymptotic UV scaling regime well above the particle 
mass and an intermediate energy scaling regime at the onset of the spinodal instability. Our 
goal is to connect the two scaling regimes and the global reconstruction of the renormalized 
trajectory. A preliminary account of our results has already been given in ref. ||. 

The possibility of the the onset of a new scaling law has been raised by the introduction 
of the "dangerously irrelevant parameters" Q. The usual classification of the operators 
in a scaling regime with respect their response to the change of the scale is based on the 
perturbation expansion and assumes the regular behavior of the system in the weak coupling 
regime. If the evolution equations develop singular behavior in the weak coupling limit 
then it may happen that a perturbatively irrelevant coupling constant experiences new, 
nonperturbative scale dependence. This mechanism can be identified in different systems 
with condensates |J, 0, M, where the singular dependence of the saddle points in the 
coupling constants provides the mechanism to turn a perturbatively irrelevant coupling 
constant into an important one. It is conjectured in this paper that the condensation might 
serve as a quite general mechanism to generate new, important coupling constants and the 
dynamical renormalization group that addresses the question of the scale dependence in 
non-equilibrium, time-dependent phenomena is well suited to the study of such a question. 

We rely in this work on the renormalization group realized by a sharp cutoff in momentum 
space, in the framework of the gradient expansion for the action. Such a realization of the 
cutoff renders the systematic gradient expansion questionable, so we restrict ourself to the 
lowest order, local-potential approximation M, |J. But it should be mentioned that the 
usual remedies of the problem, the use of a smooth cutoff or periodic Brillouin zone are 
not compatible with the loop expansion for systems with a condensate at finite-momentum 
scales. In fact, though the the loop expansion produces the action of the effective theory as 
a power series in h, the saddle point, the minimum itself is not necessarily a polynom of h. 
In order to preserve h as a small parameter, all of the stable modes have to be eliminated 
in the loop expansion before we compute the contributions of the unstable modes. Such 
a successive elimination of the degrees of freedom introduces automatically the strategy of 
the renormalization group with sharp cutoff, and renders the smooth cutoff regularization 
inconsistent with the loop expansion. 



The organization of the paper is the following. The presence of several scaling regimes is 
pointed out for the Theory of Everything and for the BCS ground state in Section 2. The role 
the Bose condensation plays in generating new scaling laws is elucidated in Section 3. The 
infinitesimal, renormalization group equation is derived for the 4 model in Section 3. The 
possibility of a new scaling regime in the vicinity of the spinodal instability is pointed out 
in Section 4. The numerical results for the solution of the renormalization group equation 
are presented in Section 5, in the case of the four dimensional 4 model. Finally, Section 6 
contains the Summary. 

II. HIDDEN COUPLING CONSTANT 

The most fundamental and at the same time the most complex appearance of the multiple 
scaling laws can be found in the renormalized trajectory of the Theory of Everything, TOE. 
The simplest procedure to describe the manner the TOE gives rise a chain of lower energy 
effective theories is the so-called matching. This method where we match two theories 
at their crossover scale can not give an account of the change of the scaling laws in a 
dynamical manner, because different sets of coupling constants are used at the two sides 
of the crossover. Instead, one should follow the original strategy of the Wilson-Kadanoff 
blocking, and put all coupling constants in the lagrangian from the very beginning that will 
later be generated by the blocking. Thus the coupling constant space of the TOE should 
contain not only the renormalizable parameters but any coupling constant we ever need in 
Physics. For example the quark-gluon vertex or a coupling constant of the Hubbard model 
of the condensed matter physics have to be considered as complicated composite operators 
in terms of the fundamental particles of the TOE. 

The renormalized trajectory, depicted schematically in Fig. [I], approaches several fixed 
points in its way towards the infrared limit. To understand this better consider the increase 
of the observational energy in the regime l-60GeV. The evolution of the coupling constants 
receive their dominant contributions from the strong interactions, from the radiative cor- 
rections of QCD, and the renormalized trajectory is in the scaling regime of QCD. Had 
our world contained the strong interactions only, the renormalized trajectory would have 
converged to the fixed point of QCD with increasing energy. But the weak interactions 
become important as we reach the characteristic energy of the Electro- Weak theory, and 
the renormalized trajectory turns away from the QCD fixed point. This happens because 
the running coupling constants of the non-renormalizable quark vertices generated by the 
exchange of the intermediate vector bosons increase with the energy. They saturate at the 
crossover between the strong and the electro-weak interactions where the guidance of the 
evolution is taken over by the exchange of the W and the Z bosons. In a similar manner, 
the fixed points of all the other, renormalizable, effective theories are approached by the 
renormalized trajectory, but the higher-energy processes always prevent convergence as the 
energy is increased, except at the last fixed point, at the TOE. In the regime of Solid State 
Physics, we can influence of the evolution of the running coupling constants by the environ- 
mental variables, such as the temperature or chemical potentials ||10|| . In this manner the 



renormalized trajectory may bifurcate and follow different path in different environments 
and finally arrive in different thermodynamical phases at the infrared fixed point. 

Facing such a complex system, the usual argument about universality appears as an 



oversimplification. In fact, at each scaling regime we classify the operator algebra of the 
model in a local manner according to the appropriate scaling laws. It may happen that 
an operator possesses different local classifications, and it is found relevant at one scaling 
regimes but becomes irrelevant at another one. The importance or unimportance of such a 
coupling constant must be decided in a scaling-regime-independent, global manner. 

A simpler example where the global behavior is important can be found in QED contain- 
ing the electron and a heavy pointlike particle with charge +Z playing the role of a nucleus. 
Suppose that certain environment variables, such as the temperature and the baryon chem- 
ical potential are chosen in such a manner that the vacuum is a solid state lattice in the 
superconducting phase. We can distinguish two asymptotic scaling regimes in this model: 

• Asymptotic UV scaling: At energies above the nucleus mass the evolution equation 
is given in terms of the minimal coupling vertices and the vacuum can be considered 
perturbative. As indicated above we do not require the existence of a fixed point for the 
identification of the asymptotic scaling laws. In this manner we can ignore the possible 
problems arising from the non-asymptotically free character of QED, and consider the 
scaling laws only up to the UV Landau pole. The relevant and the marginal operators 
are the usual renormalizable ones, and are given in the framework of the perturbation 
expansion by a power-counting argument. The size of the scaling regime is limited by 
the UV cutoff or the Landau pole. 

• Asymptotic IR scaling: For energies below the scale of eV the collective phenomena 
of the solid state lattice dominate the scaling laws. The inhomogeneity of the vacuum 
is a key element. No systematic classification of the scaling operators is known, but 
due to the massless acoustic phonons, the existence of nontrivial relevant or marginal 
operators can not be excluded. The lower edge of the scaling regime is limited by the 
loss of quantum coherence. 

The asymptotic scaling regimes can be extended by embedding the model into a higher- 
energy, more fundamental, renormalizable theory and by approaching absolute-zero temper- 
ature. These asymptotic scaling laws are important since their relevant coupling constants 
may influence the physics in a strong manner. The four different possibilities with respect 
to these scaling laws are shown in Table | and in Fig. |^ for coupling constants whose di- 
mension is removed by the cutoff, A. The electron mass is a relevant parameter of the 
QED lagrangian and remains relevant for the solid state lattice, too. The muon mass is 
relevant in the UV scaling regime but the muon-induced processes are overwhelmed by the 
electron-induced ones at low energy and m^ becomes irrelevant in the IR. The four-fermion 
interaction is, at the same time, non-renormalizable and represents the driving force to the 
BCS superconducting phase; it is marginal in the IR regime fL3| , fl4}l . Recently other indi- 



cations of the deviation from the Fermi-liquid behavior resulting from relevant or marginal 



operators of the IR regime have been found for high-T c cuprates ||16| , as well. Finally the 
six-fermion or any higher order vertex is irrelevant in the IR regime because its effect is 
reduces to the multiple application of the four-fermion vertex. The four-fermion interaction 
plays a special role: On the one hand, it is usually left out from the microscopic lagrangian 
because it is suppressed in the UV scaling regime. On the other hand, it has a key role at 
the IR regime in controlling the attraction between the electrons. 



The qualitative behavior shown in Fig. |2]c raises the following possibility. The suppres- 
sion of the irrelevant coupling constants in the UV scaling regime is used to explain the 
universal behavior of the models. It is certainly correct to expect that the renormalized 
trajectories whose initial conditions in the ultraviolet differ only in the values of the irrele- 
vant coupling constants approach each other as we move towards the infrared direction. But 
if there is another scaling regime where an operator that was irrelevant in the UV scaling 
regime turns out to be relevant, then the resulting amplification process may undo the sup- 
pression at the UV. May one find a "hidden coupling constant" in this manner that has to 
be put into the microscopic action of the UV regime but influences the dynamics only in the 
IR ? The answer to this question is nontrivial even if we can identify a nonrenormalizable 
operator that is relevant in the low-energy scaling regime. In fact, it may happen that the 
increasing value of this coupling constant during the lowering of the observational energy 
happens to be independent of its initial value at the UV cutoff. This possibility is shown 
schematically in Fig. [3|a. The UV scaling laws suppress the dependence on the intial values 
of the nonrenormalizable coupling constants. If the evolution equations have no instability 
or other nonanalytic features at low energies, then this suppressed sensitivity and univer- 
sality is observed down to zero energy. In Fig. [|b the value of the coupling constant at 
the cutoff becomes numerically negligible at intermediate energies according to universality. 
However the low-energy scaling laws amplify the numerically-small nonuniversal value. De- 
pendending on the critical exponents and the size of the scaling regimes the amplification 
of the sensitivity at the IR side may be comparable or even stronger than the suppression 
in the UV regime. This is a hidden coupling constant because it is undetectable at finite 
energies nevertheless its small value influences the IR physics in a nonuniversal manner. By 
rearranging a longer IR evolution we can in principle uncover the presence of nonuniversal 
interactions at higher energies. 

It was pointed out that a relevant operator of the IR scaling regime may change the 
ground state either by generating bound states that condense or by driving the system to 
strong couplings until the growth of the coupling constants is cut off by other quantum effects 



1 3| . Within the framework of the saddle-point expansion these two possibilities coincide, 
and the issue of the hidden coupling-constant is that the growing coupling constants might 
be slowed down in a manner that contains information about the nonuniversal interactions 
at the microscopic scale. In order to clarify this question, namely the sensitivity of the IR 
end of the renormalized trajectory on the UV initial conditions, we must obtain and solve 
the renormalization group equations for sufficiently many coupling constants, globally. 

The coexistence of several scaling regimes have already been studied in Condensed Matter 
Physics, where competing interactions are represented by the possibility of approaching 



different IR fixed points [ 15] . The difference between such systems and the TOE is that only 
one of the possible scaling laws are realized in the former case. On the contrary, the scaling 
regimes occur at different energy scales and the system visits each of them sequentially in 
the TOE. 

The renormalized trajectory of TOE sketched in Fig. [I] reflects the usual conflict between 
the "fundamental" and "applied" physics. The fundamental, microscopic parameters should 
be determined by the help of the renormalization conditions imposed at the energy scale 
Hfund that is within the scaling regime of the fundamental interactions, and is far from 
the complexity entering at lower energies. So long as the TOE is renormalizable or finite 



the perturbation expansion can be used to show that the resulting renormalized trajectory 
is independent of the choice of Hf un d, and can in a natural manner be characterized by 
the coupling constants observed at high energies. In this sense the fundamental laws of 
physics are determined by the high-energy experiments and the description of the lower- 
energy, complex systems requires "only" the capacity to apply the fundamental laws in a 
complicated situation. 

But the fallacy of this view is clear: In the absence of hidden coupling constants where 
the low-energy effective-parameters are determined in an autonomous manner, there is no 
need for the precise measurement of the high energy parameters in order to reproduce the 
low-energy physics. In fact, according to the renormalized perturbation expansion and 
universality, the bare coupling constants are characterized in a unique manner by the renor- 
malized ones defined at low energies, fi C om P i « Hfund- When hidden coupling constants are 
present then the precision of the measurement of the high-energy parameters required to 
predict a low-energy phenomenon with a reasonable accuracy might render the global de- 
termination of the renormalized trajectory illusory. In other words, the specification of the 
initial conditions with a reasonable accuracy might be unimportant for the nonlinear evo- 
lution, similarly to the nonintegrable chaotic systems. The challenge lies in comprehending 
the matching of the different scaling islands of Fig. [I]. 

III. RENORMALIZATION AND CONDENSATION 

In the perturbatively-implemented renormalization group we have to assume that the 
renormalized trajectory does not leave the vicinity of the Gaussian fixed point. If there 
is a small parameter assuring this, like in 4 — e dimensions, then the set of the relevant 
or marginal operators is the same at each scaling regime, and the global behavior of the 
renormalized trajectory contains nothing new compared to the local analysis performed at 
the individual scaling regimes. If the renormalized trajectory explores the nonperturbative 
regions then the global analysis is very interesting, but we loose the general, perturbative 
characterization of the flow. A window of opportunity for the analytical studies opens if the 
nonperturbative features of the evolution can be reproduced by the only known, systematic, 
analytical, nonperturbative method, the semiclassical expansion. 

A nontrivial saddle point corresponds to a coherent state formed by the condensate of 
bosons. The importance of this condensate depends on the multiplicity of the coherent 
state(s). We shall demonstrate this from three different points of view. 

• Euclidean Quantum Field Theory for the vacuum. In the semiclassical solution to 
the vacuum with ferromagnetic condensate only the lowest lying excitation level with 
vanishing momentum is populated macroscopically. [] When the system is placed in a 
finite geometry then the lowest lying state might become inhomogeneous; nevertheless 
the particles still condense in a single mode. The impact of the condensate is stronger 
if the particles form coherent states in a large number of modes. This is the case for 



x The degeneracy of the vacuum with spontaneously broken symmetries does not change the picture 
since there is only one vacuum in each "world" containing the states connected by local operators. 



solitons or instantons, the localized saddle points with high entropy. Another, less 
studied example is given by theories with higher-order derivative-terms in the action 
M. When their coupling constant is properly chosen then particles with nonvanishing 
momentum may condense and generate an inhomogeneous vacuum. Since the density 
of modes with momentum p is proportional to p d_1 in d dimensions there are more 
particle modes participating in the condensation at finite p than at p — 0. 

Real time dynamics. The condensation is triggered by the wrong sign of the forces 
trying to restore the equilibrium position. There are always several unstable modes 
which compete in the condensation when the system starts with the naive vacuum 



njj. Depending on the initial conditions different slower collective modes might be 
formed that distribute the energy and create time dependent condensates in several 
other particle modes. After sufficiently long times, the energy of the excited modes get 
diffused by the friction terms over infinitely many modes, without creating macroscopic 
population, and only the particles in the lowest lying mode remain condensed. The 
large- amplitude, transient excitations above the true vacuum are described by a large 
number of modes with condensate. This is the realm of nucleation and the spinodal 



decomposition [12 



The former refers to the metastablility of the vacuum which decays by large am- 
plitude fluctuations. The simplest is to assume that this takes place by the spon- 
taneous formation of spherical droplets of the stable vacuum whose free energy is 
F = 4nRa — 4nR 3 AF/3, where a, AF and R stand for the surface tension, the 
free energy density difference betwen the metastable and the stable vacuum and the 
droplet radius, respectively. These droplets extend over the whole volume if they are 
sufficiently large, R > R cr = a/AF. The spinodal decomposition is observed when 
there is no more finite threshold for the instability and the infinitesimal fluctuations 
are enough to trigger the decay of the homogeneous vacuum. This can be recognized 
in the framework of the static description by the apparence of the negative eigenval- 
ues for the Euclidean propagator. The fast increase of the corresponding elementary 
excitation amplitudes drives an inhomogeneous separation of two stable values of the 
local dynamical variable. 

Returning to the long-time, low-energy excitations, they experience a single condensed 
mode; and the static nonperturbative condensate is simpler to describe. 

Saddle points of the renormalization group. When the renormalization group is im- 
plemented then we eliminate the modes in descending order in the energy, and we 
may encounter nontrivial saddle points in the way. If we are interested in the effective 
theory for the low-energy fluctuations around the true vacuum in systems without a lo- 
calized saddle point, then the condensate occurs only at the last mode, at the infrared 
fixed point. When the effective theory is sought for large amplitude excitations, then 
the eliminations of the modes in the presence of such a background field may induce 
saddle points for the blocking procedure earlier, at finite scales. The nonperturbative 
contribution of these saddle points may modify the direction of the renormalization 
group flow in a substantial manner [ 17 1 . As an example of this mechanism, the emer- 



gence of an irrelevant coupling constant in the dynamics of domain walls has already 
been noted in the two dimensional nonlinear 0(2) model Jl8| . 



In short, the saddle points of the blocking procedure occur with higher multiplicity at 
finite energies, and may modify the scaling laws in a more substantial manner than the saddle 
points of the static system in the true vacuum. Such effects could be seen by the proper 
implementation of the dynamical renormalization group for large amplitude fluctuations. It 
is worthwhile noting that the spinodal instability region starts just at the vacuum in the 
presence of the Goldstone modes, and the saddle point contributions might influence the 
small quantum fluctuations around the true vacuum, as well [T3 . 



IV. RENORMALIZATION GROUP EQUATION 

The study of the different scaling regimes requires the handling of a large number of 
coupling constants which can be achieved by converting the evolution equations referring to 
the individual coupling constants into a differential equation for the generating function for 
the coupling constants, i.e. the blocked action. We shall work with the four- dimensional 
4 theory in Euclidean space-time. The cutoff k in the momentum space will be changed 
infinitesimally, k — > k — Ak, which generates a new small parameter, Ak/M, where M is 
a quantity of dimension of mass made up by the dimensional parameters of the model and 
the cutoff, k. This small parameter will be used to show that the renormalization group 
equation derived in the one-loop approximation becomes exact as Ak/M — > 0, because the 
higher loop contributions are suppressed by Ak/M. 

There are two limitations to bear in mind in turning this scheme into a feasible algorithm. 
The first is that even thought there is a new small parameter to suppress the higher loop 
contributions to the evolution equations the argument presupposes the applicability of the 
loop expansion. Therefore the use of the resulting "exact" equation is questionable in the 
strong coupling region, beyond the validity of the loop expansion. The other limitation 
comes from the parametrization of the effective-action functional, S k [(f)], with the cutoff k. 
It is usually done my relying on the gradient expansion, 

1 



S k [(f>] = Jd d x 



ZMx^d^x)) 2 + U k (c/>(x)) 



0{d% (1) 



The assumption that terms with higher order derivatives are less important is equivalent 
with the belief that the action is a local functional. This is reasonable at high energies 
because all of the relevant operators of the short distance scaling laws should be local. But 
there is no reason to exclude operators which are nonlocal at the scale of the cutoff at 
low-energy, scaling regimes, and these terms of the action may create nonlocal effects. By 
holding to the assumption of locality we shall set 

Z k {<j>) = 1 (2) 

in this work. 

We eliminate the modes with momentum k — Ak < p < k and find the blocked action 

^Sk-AkW = f D[0 / ] e -i s *^+^'] (3) 



e ^ 

The Fourier transforms of the variable <fi{p) and 4>'{p) are nonvanishing for p < k — Ak and 
k — Ak < p < k, respectively. The functional integration is carried out in the framework of 
the loop expansion, 

8 



S k -Ak[<P] = S k [(f> + <f>' ] + hr'ln5 2 S + 0(h 2 ), (4) 

where tr' denotes the summation within the shell [k — Ak,k}. The positive semi-definite 
operator 

p S(x v) - gg#±fj (5) 



and the saddle point, <f/ 



o- 



ssS + P] = (6) 

W( X ) \4>>{x)=4>' {x) 

are computed by keeping the background field, (f>(x), fixed. Eq. (|) reduces to the Wegner- 
Houghton equation || when the saddle point is trivial, <p' — 0. 

Observe that the loop integration is made in a region of volume Vt ( ik d ~ 1 Ak where 

2 d ir d/2 
Qd = f(d/2) " (7) 

The contributions 0(h n ) are given in terms of n-fold loop integrals. So long as the integrands 
are bounded in the domain of integration the contributions 0(h n ) contain the multiplicative 
factor (Ak/M) n . The integrands are the products the propagator, the inverse of (||) evalu- 
ated within the subspace k — Ak < p < k. If the restoring force acting on the fluctuations is 
nonvanishing then the propagator is bounded and the higher loop contributions drop from 
(||) when Ak/M — > 0. The formal argument showing the suppression of the higher loop 
contribution to the renormalization group equation is presented in the Appendix. 

The only quantity needed in the local potential approximation, (fO)-(0), is the potential 
£/fc(<&) so it is natural to choose for its determination a homogeneous background field, 
4>(x) = $, for which the kinetic energy is vanishing. An additional bonus of this choice is 
that the saddle point is nontrivial, <p' ^ 0, just in the spinodal unstable phase. Outside of 
the spinodal unstable phase one can set <j)' Q = what yields 



e 



5fc[*+^o]-|tr' ln5 2 S+0(h 2 ) 



= e -L d U h m-^f *faH[j?+Ui!(*j\+0{h*) ^ , g v 

where the plane wave elementary excitations of the homogeneous vacuum were used to 
compute the trace in the third line and J d A p = fk-Ak<p<k ^V- in this manner one arrives 
at the finite difference equation 



U k (<S>)-U k - Ak (<S>) = --J-^- d \n{p 2 + Um 



(2tt) c 
where 



1 + (^)j 



(9) 



<*<•) - ^ m 



By taking the limit Ak ->Owe find the differential equation 

*>'<•> = -»$&'■[*' + '*<•>]■ i 11 ' 

by the help of the polar coordinates in the momentum space. 

The evolution equation, ([11]), is actually the one-loop resummation of the perturbation 
expansion. This can be seen by expanding in the non-Gaussian pieces of the potential, 



k—U k (<£>) = -h 



d 



dk " v ' 2(2vr) 



ln(A: 2 + ml) + £ 



\n+l 



'U»(®)-ml 



,,-i n \ k 2 + m\ 



(12) 



where m 2 (k) = U{!(Q). The evolution equation for the effective vertices at zero momentum, 
i.e. the leading order of the gradient expansion is now the explicit sum of the one-loop 
graphs defined by the free propagator G~ 1 (k) = k 2 + m|. If we ignore the fc-dependence in 
the right hand side then Uo(<0>) is the usual one-loop effective potential obtained in the bare 
or the renormalized perturbation expansion for 

m\k)-,m\ UZ{*)-ml-+?£&, (13) 

or 

m\k)^ml &»(<&) - m *-^#», (14) 

respectively. The evolution generated by the renormalization group provides a partial re- 
summation of the perturbation expansion by piling up the effects of the modes which have 
already been eliminated in the /c-dependence of the running coupling constants, £/&(<&), in our 
case. In the usual renormalization group method, which is based on the renormalized pertur- 
bation expansion, only the renormalizable coupling constants, d 2 Uk{0)/d<& 2 , d 4 Uk(0) / d& A , 
pile up the effects of the elimination. In solving the differential equation ( |i"2"|) the nonrenor- 
malizable coupling constants, i.e. the higher order derivatives are evolved, as well. Such an 
extension of the usual scheme is necessary to find the eventual hidden coupling constants. 

It is easy to perform the direct numerical integration of (|ll|) with an initial condition 
U\(Q) given at k = A towards the infrared direction after sufficient attention is paid to 
the instabilities arising from the finite $ and k resolution. However our goal is to arrive 
at a description of the modification of the scaling laws which requires that the restoring 
force for the fluctuations, the argument of the logarithm function in (|TTD , be vanishing at 
k = fco($) 7^ 0. The singular behavior of the differential equation made the numerical 
quadratures unstable and forced us to resort to methods other than the discretization of the 
$-space. Instead of discretizing the variable $ we truncated the Taylor expansion of Uk(4>) 
at finite order. 

To this end we introduce the coupling constants defined at $ = $ by 

The result is 

10 






The /3-functions are defined as 



Pn = k-g n (k) = Q^k m U h (*) 



<j>=<V 



(16) 



(17) 



where the analyticity of the potential in k and $ was assumed. They are obtained by taking 
the successive derivatives of (|TT|), 



"- - - s t5^< G " 



Gn 



+2) 



where 



and 



V 



G, 



gn 



fjn 



k 2 + g 2 ' 



<9$ n 



In 



k 2 + wim 



(18) 



(19) 



(20) 



is a polynom of order n/2 in the variables Gj, j — 2, ■ • • , n + 2. 

In order to exploit the simplification offered by the symmetry £/&(—$) = [/&($) we set 
$o = 0, which cancels the odd vertices and yields for the even ones 



V 2 = Gi, 



i- 



(21) 



rQ — Cxg — 15GgGr4 + o0Gr4, 



One can verify diagrammatically that the system of equations (pO|)-(^l]) is a compact rewrit- 
ing of the one-loop contribution to the beta functions of the effective vertices. 

It is useful to obtain the evolution equations for the coupling constants whose dimension 
is removed by the cutoff, 



Pn = k-^rg n (k) 



d n - d ~ - 
d$ n dk k 



(22) 



4U 



d 



~ 9n ~ n 2^r Vn{ ^ 



1 G n ^ 



where 



and 



k = AJfe, U = k d U, $ = k-2- 1 ®, g n = k n{1 -i )+d ~g n , 



G, 



9n 

1 + &' 



(23) 



(24) 



11 



V. SCALING REGIMES 

The 4 model has two scaling regimes in the symmetrical phase, separated by a crossover 
at the mass, k 2 r = g2{k 2 T ). For k » k cr 



Gn ~Y 2 



\k 2 



G n = g n [l + 0{g 2 )\, (25) 



and the evolution of the coupling constants is those of a massless theory with the only scale 
of the cutoff. Below the crossover scale we have 




G n = ^[l + 0(~g^)], (26) 

#2 



n _9n 

92 _ 

and the evolution comes to a halt due to the factor k d in right hand side of fllT]). This is as 
expected since a theory with mass gap has the quadratic mass term as relevant operator and 
all non-Gaussian coupling constants are irrelevant. The /3-functions are always dominated 
by the highest order coupling constant that is linear in the coupling constant in question. 
In fact, this term contains the lowest power of 1/k 2 at high energies. At the IR side this 
term is the most important because the others contain the product of more non-Gaussian 
coupling constants which are supposed to be small. 

Qualitatively new scaling proprieties can be established if the /3-functions are dominated 
by other terms or receive comparable contributions from different terms. This is certainly 
the case when the restoring force, 

D(k) = k 2 + g 2 (k), D(k) = l + g 2 (k), (27) 

is vanishing. Suppose that a theory specified at the cutoff by U\(&) is in the symmetrical 
phase, i.e. g 2 (0) > for the choice $ — 0. The mass gap decorrelates the field variables 
which are well separated in the space-time, and the central limit theorem asserts that £4($) 
approaches a quadratic form and g 2 (k) is monotonically increasing as k — > 0. Perturbation 
expansion gives a more detailed picture: g 2 (k) is a monotonically decreasing {fi 2 < 0) and 
D(k) > is a monotonically increasing function for all values of k. In order to find the 
new scaling laws we need either a massless or symmetry broken-theory, where D(k) reaches 
zero at k = or at k = k ^ 0, respectively. In the former case the Coleman- Weinberg 
mechanism [BOj generates a different vacuum where the theory manages to develop a mass 
gap, thereby preserving the usual IR scaling laws. In the latter case the vanishing of the 
inverse propagator, D(k), indicates an instability of the homogeneous vacuum in the effective 
theory with the cutoff k <k$ with respect to small fluctuations, the presence of the spinodal 
instability. Thus new scaling laws might be found as the precursor of the spinodal phase 
separation. We should bear in mind that within the spinodal-unstable region the saddle 
point is nontrivial, <j)' ^ 0, and we need a different renormalization group equation. 



VI. RENORMALIZED TRAJECTORY 

We want to follow the renormalized trajectory of the coupling constants 
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92n{k) ~ d^ |* 0= o' (28) 



in d = 4, i.e. we seek for the solutions of the set of coupled equation 

d_ 

'dk 



d 

k—g n (k)=l3 n (g u ---,g n+2 ), (29) 



where 

,4 ^n 

'2(270* 0$ n 

with the initial condition 



Pn = -n^^ln\k* + U'a*)l, (30) 



$o=0 



U A ($) = U B ($) (31) 



in the stable region, 



*?>%{*) = -Ui, m {*). (32) 

We set h = 1 in the numerical work. We truncated the potential 

W) = E|^ 2n , (33) 

and solved the resulting equations numerically. Special care is needed in the vicinity of the 
spinodal line, k ~ &o(<&), where the /3-functions are the sum of large numbers with different 
sign. We used 3-nd and 4-th order Runge-Kutta method with dynamically determined 
value of Ak. Quadruple precision numbers were used when necessary to make sure that 
the roundoff errors in the /3-functions were less than 10 -8 times the actual value of the 
/3-functions. 

The behavior of the restoring force for the fluctuations, D(k), was found in qualitative 
agreement with the perturbation expansion, two typical cases are shown in Fig. f|. Such 
an agreement is expected in the UV scaling regime where the kinetic energy dominates the 
action. Our interest will be to see the detailed behavior of the flow in the vicinity of the 
critical curve, k > k Q (Q). The curves /co($) obtained from the tree-level solution and the 
numerical integration of the renormalization group equations are shown in Fig. [5[ The 
quantum fluctuations help the disorder and drive the saddle points to zero at the weakly 
unstable regime of the tree-level scaling relations. As a result the curve obtained by solving 
the renormalization group equations is inside of the spinodal-unstable region of the tree 
level solution, except for the roundoff effect at small k. It is remarkable that the radiative 
corrections are rather small. 

The numerical integration of the renormalization group equations produces oscillation 
for the coupling constants g n (k) for n > 4, with increasing amplitude as we approach the 
instability. We found two qualitatively different behaviors as far as the vicinity of the 
unstable line is concerned. 

Focusing. When the potential is truncated up to 2N = 20 then gjjk) stays positive and 
approaches zero. The values of the higher order coupling constants drop significantly after 
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several, large-amplitude oscillations and approach zero. This indicates that the blocking 
transformations has an attractive fixed point, 

The coupling constants g n (k) for n > 4 produce large fluctuations, but after a while start 
to all fall and the quadratic potential, (0) is approached. It is found that g±{k) — > + and 
the higher order coupling constants drop after undergoing large amplitude oscillations as 
k — > feo- The evolution of In \g2o{k)\ is depicted in Fig. ^. We find a cusp where the sign of 
92o{k) changes with finite altitude due to the finite resolution of the k values. Note that the 
negative coefficient of <fi 2 causes no problem with the stability of the vacuum because the 
potential £4($) of the effective theory recovers the perturbative form for large values of the 
field, far away from the spinodal-unstable region. 

The approach to the Gaussian potential can be made plausible by inspecting V±. By 
assuming that the coupling constants remain finite at the critical line we have ^(^(O)) = 0, 
since (3^ diverges otherwise. Once g^ = is accepted, the vanishing of the higher order 
coupling constants is plausible. The only finite parameter, #2(^0(0)), is fixed by the condition 
D(ko(0)) = 0. The existence of a single fixed point, that all finite set of coupling constants 
runs into Q54|), can be called focusing. The potential turns out to be ( |54"D in the whole 
unstable region that appears as a "fixed region" ||19|| . 

Divergence. This fixed point turned out to be an artifact of the truncation of the poten- 



tial, a feature what has already been noted in other cases, as well, |2l|. When the truncation 
is made beyond 2N = 20 then the accumulation of the contributions of the higher order 
vertices in the /3-functions make the term G$ more important in (3^, which in turn helps g± 
to decrease faster with k in the approach of the critical line. Once the sign of #4, flips the 
further decrease is not limited by zero and #4 quickly approaches — 00. It is not then so 
surprising that all of the other coupling constants start to diverge at the same time. The 
typical flow is depicted in Fig. [7] for 2N = 22, the further increase of N makes no further 
qualitative changes in the flow. The coupling constant undergoes oscillations with increasing 
amplitude as g^{k) — > + and start to diverge as g&{k) flips sign. The loop expansion natu- 
rally ceases to be applicable in the vicinity of the critical line, and all we can say is that the 
modes with momentum slightly above /c ($) appear strongly coupled, and our solution is no 
longer reliable. Though there is a marked difference in the behavior of the renormalization 
group flow for 2N < 22 and 2N > 22 but we should not forget that this difference shows 
up after a strong coupling regime where the high order non-Gaussian coupling constants 
develop extremely large values. So it is not clear if the difference between two two cases in 
the vicinity of the unstable line is indeed so large. 

Universality. The increase of the coupling constants at k « k (0) indicates the existence 
of new relevant operator(s) in this scaling regime. One suspects that this operator is nonlocal 
since the value of the cutoff is finite. Can this operator modify the usual universality 
argument || ? According to universality, the introduction of the irrelevant operators at the 
cutoff modify only the scale parameter of the theory. The dimensionless quantities, such 
as the /3-functions, are supposed to be independent of the value of the irrelevant coupling 
constants at the cutoff. To verify this scenario we computed 

%w (35) 

14 



numerically. The result, plotted in Fig. |8] shows clearly the coexistence of two different 
scaling regimes, the UV one where this quantity is suppressed and the precursor of the 
spinodal instability where we find an increasing value. 

The comparison of Figs. [7] and § contains an important lesson confirmed by Fig. |9| where 
(|35|) is plotted as the function of the appropriate coupling constant. Namely, the violation of 
universality, the increase of ( |35|) takes already place when the coupling constants are weak. 
It is reasonable to assume the ansatz 

998(1) ^ k (0) ) ayh 

where the term /c 2 (0) that is proportional to A -2 represents the suppression of the UV scaling 
laws, and the amplification effect of the instability is manifest in the behavior of F(x), that 
is supposed to diverge after large oscillations, as x —* 0. For any finite value of the cutoff 
we can find a value of k sufficiently close to the instability where (^) is unity, 

and the cutoff effects of the nonrenormalizable coupling constant become visible because 
the suppression of the UV scaling regime is compensated for by the amplification of the 
instability region. 

If the environment of the system, represented by the insertion of the constraint 

5 U - - f d 4 x<f){x)) (38) 

where V is the four volume into the path integration is chosen in such a manner that the 
spinodal instability occurs, i.e. $ is within the spinodal unstable region, then the amplifi- 
cation makes the effective coupling strength nonuniversal around the instability. In other 
words, the coupling constants of the effective theory for the system subject to this constraint 
can pick up the the values of certain nonrenormalizable coupling constants at high energies, 
and allow us to investigate the nonuniversal interactions at high energy if the cutoff of the 
effective theory is brought close to k — fco($)- By this method one could in principle increase 
the energy regime we can access experimentally, and may get closer to the " last" important 
scale, the onset of the asymptotic UV scaling of the TOE. This mechanism can be called 
a "renormalization group microscope" since the amplification offered by the instabilities is 
similar to the usual microscope, except that it is achieved by the renormalization group flow 
in the space of the coupling constants. 

According to the numerical results the strength of the singularity and the value of ( pop 
approach zero at the critical line, as $ is increased towards the edge of the unstable region 
k = ko(<&). This is the result of the factor k d in the renormalization group equations, the 
decreasing entropy of the modes with weak restoring force as k — > 0. 

Renormalized perturbation expansion. In Statistical Physics one usually follows the evo- 
lution of the bare coupling constants, g n (k), as the functions of the cutoff, 

k^g n (k)=/3i B \{g}), (39) 
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where the explicit dependence on k drops from the bare /3-functions in the UV region. We 
studied the evolution of these coupling constants in our work, too. 

In Particle Physics one introduces the renormalized running coupling constants, A n (/x), 
n = l,---,n r , where n r is the number of renormalizable parameters in the theory. The 
irrelevant coupling constants are neglected because the cutoff is sent sufficiently far from 
the observational energy, \x. The running coupling constants are defined by the help of the 
scattering amplitudes or Green functions and their scale dependence is described by the 
renormalization group equations, 



(40) 



involving the renormalized /^-functions. 

So long as the perturbation expansion is applicable we can establish a one-to-one mapping 
between these schemes, 



9n(k) = X n (k) + 



(41) 



where the terms omitted are higher orders in the supposedly small coupling constants or 
are 0(k/A). In this manner the qualitative features of the renormalization group flow agree 
in the two schemes, i.e. the running coupling constant becomes the bare one when the 
scale reaches the cutoff. How can we reconcile the importance of certain nonrenormalizable 
coupling constants at low energy with the evolution of the running coupling constants, A n (/x), 
which are introduced to keep track of the renormalizable coupling constants only ? 
The independence of the theory from the observational scale, 



/^r B ({A B },A) = 0, 



(42) 



is reached by the readjustment of the coupling constants, A„(/i). According to the multi- 
plicative renormalization scheme the number of adjustable parameters is just n r + 1. In 
particular, one chooses n r + 1 independent observables, r( m ) ; rn = 1, • • • , n r + 1, 



rt } ({A B },A) = Zh (W)},^) 4 m) ({A(^)},/x), 



and imposes 



d 
= /i— 
d\i 

d 



^(W)},£)ri m) ({A(^)},^ 

d 



^ + ^ R) (i9})-^- n + tmi {R \{g})\ rP({\(f,)},v). 



(43) 

(44) 
(45) 



The renormalized renormalization group functions, 7^ and /3^({A}) can be found by 
inverting 



= n 



({\(n)},f\rP({\(n)},n) 



djj, 



C7 (R) (M) V { r\{K»)}^) 



(46) 
(47) 
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for the 7^ and the @( R ' functions. The perturbative renormalization assures that these 
functions are well defined, i.e. are independent of the choice of the observables {T}. 

Observe that there are two steps in this procedure which prevent the detection of a 
hidden coupling constant. One is that we have already committed ourselves to use as 
many coupling constants to compensate the modification of //, as many renormalizable 
operators are in the system. The second is that the cutoff is formally removed from the 
renormalized perturbation expansion. The hidden coupling constant represents a coupling 
between the ultraviolet and the infrared modes that can be seen by keeping the UV and IR 
cutoffs finite. When the UV cutoff is removed then the contributions 0(ji/A) are ignored. 
The renormalization group used in this work resums these contributions which yield the 
new nontrivial scaling laws. Thus the observables obtained by the improved renormalized 
perturbation expansion using the renormalization group for n r coupling constants will, by 
construction, never show any indication of the eventual, hidden parameters. What we have 
found is that our partial resummation of the perturbation expansion, which keeps track of 
nonrenormalizable operators as well, indicates the presence of hidden coupling constants and 
suggests that the low energy dynamics is parametrized by more than n r coupling constants. 



Consequently more than n r + 1 observables must be used in fl47|) to obtain the evolution of 
the parameters that can compensate the change of the observational scale. 

In order to find out the number of real parameters, we need a controllable method to 
study the low-energy scaling behavior. Due to the limitation of the gradient expansion based 
on local operators we can not at the present stage clarify this point. 

VII. SUMMARY 

The renormalization group is traditionally used to follow the scale behavior in the vicinity 
of a fixed point of the blocking transformation. We showed in the case of the 4 model that, 
after paying the price of following the mixing of a large number of operators during the 
blocking, the investigation of the manner by which the different scaling regimes give rise to 
each other is feasible. 

A new finite energy scaling regime of the 4 model with spontaneously broken symmetry 
is generated by the spinodal instability. It was found, by the numerical integration of the 
Wegner-Houghton equation in the local-potential approximation, that the spinodal instabil- 
ity generates new relevant operators, and may undo the suppression of the nonrenormalizable 
operators at the UV scaling regime. A nonrenormalizable operator gives rise to a hidden 
coupling constant if the operator in question is relevant in the low-energy scaling regime, 
and the initial, high energy value of its coupling constant influences the low energy physics. 
This raises the possibility of the eventual use of this instability as a renormalization group 
microscope to detect the nonuniversal physics at high-energy by going sufficiently close to 
the unstable region. 

The instability studied in this work appears in the Euclidean effective-theory at finite- 
momentum scales. It is conjectured that the simplest manner to observe the effects of 
the instability is in the framework of the dynamical renormalization group, applied for 
large amplitude fluctuations. The saddle-point structure of the effective theory, for the real 
time dependence, is needed to make more definite proposal, concerning the whereabouts 
of nonuniversal phenomena. The investigation of the renormalization group equation with 
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nontrivial condensate within the spinodal phase separation, and the search for the effects of 



the instability in the vacuum with Goldstone modes, is in progress [19 
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FIGURES 




FIG. 1. The renormalized trajectory of the Theory of Everything (TOE). It passes by the fixed 
points of the Grand Unified Models (GUT), the unified Electro- Weak theory (EW), the strong 
interactions (QCD), the electromagnetic interactions (QED) certain fixed points of the Solid State 
and Condensed Matter Physics (CM) and finally approaches the ultimate IR fixed point. The 
trajectory may be influenced by the environment and reach different thermodynamical phases in 
the IR regime. 
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(d) 
FIG. 2. The qualitative dependence of the running coupling constant of Table 1 as the function 

of the cutoff, a = 2tt/A. The asymptotic UV and IR scaling regimes are shown. The coupling 

constant is supposed to be constant in between for simplicity. 
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(b) 
FIG. 3. Two possible behaviors of a nonrenormalizable coupling constant what is relevant at 

the low energy scaling regime, case (c) in Table 1. Each plot shows two curves what belong to 

initial conditions at the cutoff what differ in the value of the irrelevant coupling constant only. 
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FIG. 4. The evolution of the inverse propagator, D(k) obtained by 2N = 22. (a) Symmetrical 
phase, 52(1) = 0.1, 54(1) = 0.01, (b) Symmetry broken phase, §2(1) = —0.1, 54(1) = 0.01 In both 
cases <72n(l) = 0, for n > 2. 
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FIG. 5. The line of singularity of the renormalized trajectory on the plane ($, k) for 52(f) = — •!> 
54(f) = 0.2. The diamond and the cross show the tree level and the renormalization group results. 
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FIG. 6. The evolution of the coupling constant g2o(k) for 2N = 20, with the initial conditions 
1) = -0.1, <j 4 (l) = 0.01 and g 2n {l) = 0.0 for n = 3, • • • , 10. 
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FIG. 7. The evolution of: (a) In | ^4 ( fc ) | , (b) In |^6(fc)| and (c) In |^22 (^) I at 2N = 22, with the 
initial conditions 52(1) = —0.1, 54(1) = 0.01 and #2n(l) = 0.0 for n = 3, ■ ■ ■ , 11. The coupling 
constants oscillate for n > 4 with increasing amplitude and changing sign so In |g„(A;)| is plotted. 



30 




31 




32 



160 



140 



120 



100 




FIG. 8. The evolution of In |^{§| with 2N = 22. (a): n = 4, (b): n = 6, (c): n = 22. 
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FIG. 9. hx\dPn(k)/dge(l)\ plotted against g n {k) for (a): n = 4, (b): n = 6, (c): n = 22, 
0.317 < k < 1. As the scale parameter A; decreases we move towards left or right for (a) or (c), 
respectively. The derivative of the beta function shows the onset of the low energy, nonuniversal 
scaling already when the coupling constants are still weak enough to rely on the loop expansion. 
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TABLE I. The four classes of the coupling constants in QED. 
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APPENDIX A: WEGNER-HOUGHTON EQUATIONS 

The heuristic derivation of the renormalization group equation (|TT|) indicates that the 
higher loop contributions to the equation are suppressed. This is not obvious from the 
derivation presented above because by placing the system into a finite quantization box the 
spectrum of the momentum becomes discrete and we may eliminate the modes ony-by-one. 
What is the small parameter in this case? 

The expansion of the action around the constant background $ in powers of the Fourier 
components of fluctuations 0' is 



dS k 



^+^'] = 54$]+e^i^,e^; 



d 2 S k 



+ 



(Al) 



p.q 



We have 



OS, 



d 2 S k 



d n S k 



, = u£\<b)L d 5{p) 
v 

U { k 2 \<$>)+p 2 ]L d 5( Pl +p 2 ) 
ut ) ^)L d 5( Pl + ...+p 2 ) 



(A2) 



y pi ...^YPn 



where the subscript stands for the derivatives with respect to $. Since k — Ak < \p\ < k, 
the first derivative of the action does not contribute, 



S k [<5> 



S k [$] 



E^-p u^w+p 



+ 



(A3) 



The minimum value of Ak being 2ir/L where L is the length of the quantization box, the 
number of modes to eliminate in the shell k — Ak <\p\ < k is 



d-l 



*< = lMW = to[to) LAk 



(A4) 



where fid is the solid angle in dimension d. The integration over degrees of freedom <$' 
will be done after the expansion of the exponential around the free action. The only terms 
contributing in the integration are those for which the Fourier components of 0' are combined 
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where V[<p'] = Y[k-Ak<\ P \<kd'R-e((j)' p )d2m((f)' p ). The Gaussian integrations lead to 

\M d /2 



-V(^-A*(*H»(*)) 



hn 



1 



We now introduce the small variable, 



k 2 + ui 2) {$)J 

h [/f(<&) 



(A6) 



^[P + ^ 2) (0)] 2 



M d 



M d 



i + 



r (2n) 



-,n-l 



f/r J (0)i wo 



[p + f/f(0)] 2 (^ 



cHn-1 



« 1 



(A7) 



what is proportional with the ratio of the number of degrees of freedom eliminated in a 
blocking to those left in the effective theory. We must keep in mind that first we have to 
take the thermodynamical limit, L —>■ oo and after Ak — > to make sure that the higher 
loop contributions are small. In principle Ak has a lower bound 2ir/L but we can imagine 
that we make an interpolation on the renormalized trajectory so that the value of Ak we use 
to derive the renormalization group equation is as small as we wish. In this manner Ak and 
L are independent and the small parameter is Ak/M. If we want to keep the lower bound 
for Ak nonvanishing then we must assume that the derivatives of the potential are small 
enough according to ( JA~7j ) and therefore the existence of another small parameters hidden 
in the potential. 

By assuming that ( |A7| ) is valid and taking the logarithm of both sides of (|A6|) we obtain 
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what yields finally to the Wegner-Houghton equation 
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where the denominator in the logarithm function was inserted to cancel the potential at 
$ = 0. 
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